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The objective of this study is to identify and characterize the 
genetic variants related to the glomerular filtration rate (GFR) 
linkage on 2q37. Of the positional candidate genes, we selected 
IRSl and resequenced its 2-kb promoter region and exons for 
sequence variants in 32 subjects. A total of 11 single nucleotide 
polymorphisms (SNPs) were identified. To comprehensively 
cover the 59-kb-long intron-1, eight additional tagging SNPs were 
selected from the HapMap. All the 19 SNPs were genotyped by 
TaqMan Assay in the entire data set (N = 670; 39 families). Asso- 
ciation analyses between the SNPs and GFR and type 2 diabetes- 
related traits were performed using the measured genotype ap- 
proach. Of the SNPs examined for association, only the Gly(972) 
Arg variant of IRSl exhibited a significant association with GFR 
(P = 0.0006) and serum triglycerides levels (P = 0.003), after 
accounting for trait-specific covariate effects. Carriers of 
Arg972 had significantly decreased GFR values. Gly(972)Arg con- 
tributed to 26% of the linkage signal on 2q. Expression of IRSl 
mutant Arg972 in human mesangial cells significantly reduced the 
insulin-stimulated phosphorylation of IRSl and Akt kinase. Taken 
together, the data provide the first evidence that genetic variation 
in IRSl may influence variation in GFR probably through im- 
paired insulin receptor signaling. Diabetes 61:2385-2393, 2012 




Diabetic nephropathy (DN) is a serious micro- 
vascular complication of both type 1 and type 2 
diabetes (T2DM). Elevated urinary albumin ex- 
cretion and decreased glomerular filtration rate 
(GFR) are risk factors for cardiovascular morbidity and 
mortahty and progression to end-stage renal failure in 
individuals with diabetes. End-stage renal disease second- 
ary to DN accounts for —50% of cases in most developed 
countries. The pathogenesis of DN is multifactorial. Fa- 
milial clustering of DN and related traits such as albumin- 
uria and GFR indicate that genes play a major role in its 
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susceptibility (1). Despite the recent efforts to identify and 
characterize DN susceptibility genes, the genes involved in 
susceptibility to the development and or progression of DN 
and in particular GFR have yet to be identified. 

Equations that estimate GFR facilitate the diagnosis, 
evaluation, and management of chronic kidney disease. 
GFR varies between individuals and is influenced by ge- 
netic and environmental factors and their interactions. Both 
genome-wide linkage and association analytical appro- 
aches have been used to localize the susceptibility genes/ 
variants that influence variation in GFR (2-11). However, 
the functional signiflcance of the flnding from such studies 
is yet to be explored. To localize genes contributing to 
variation in GFR, we recently performed a genome-wide 
linkage screen using GFR data (N = 453) from the Mexican 
American participants of the San Antonio Family Diabetes/ 
Gallbladder Study (SAFDGS) and identifled the strongest 
linkage of GFR to occur on chromosome 2q35-37 near the 
markers D2S1363-D2S427 (LOD = 3.8; LODc = 3.3) after 
incorporating genotype by environment (T2DM) interac- 
tion effects (5). Identiflcation of linkage of GFR on 2q is 
detailed in the Supplementary Data. We hypothesize that 
there is a gene(s) located on chromosome 2q35-37 that 
may be related to the linkage signal identified. 

Of the positional candidate genes in the critical interval, 
IRSl gene is located 200 kb telomeric to the marker 
D2S1363. When phosphorylated, IRSl acts as a docking 
protein for multiple Src homology-2 (SH2)-containing pro- 
teins such as phosphatidylinositol-3-kinase (PI3K), growth 
factor receptor-binding protein-2, and tyrosine phospha- 
tase SHP2. These proteins in turn further activate down- 
stream effectors that mediate the biological effects of 
insulin in several tissues including kidney (12,13). Mice 
deficient in IRSl by targeted disruption of its gene display 
insulin resistance and impaired glucose tolerance (14), and 
mice heterozygous for defects in genes for both IRSl and 
the insulin receptor (15) or IRSl and glucokinase (16) 
develop overt diabetes. Genetic variation in IRSl, partic- 
ularly Gly(972)Arg (rsl801278) variant, plays a contribu- 
tory role in insulin resistance and diabetes susceptibility 
(17) probably through impaired insulin receptor signaling 
(18-20). Genetic variants (rs2943641; rs2943650) located 
500 kb upstream of IRSl were recently reported to be 
associated with insulin resistance, hyperinsulinemia, dys- 
lipidemia, and T2DM (21,22). 

Given the functional significance of IRSl together with 
our finding that it is a positional candidate gene in the 
GFR linkage region on chromosome 2q, we determined 

DIABETES, VOL. 61, SEPTEMBER 2012 2385 



IRS1 GENE AND GFR 



whether genetic variation in the IRSl gene contributes to 
variation in GFR in the Mexican American participants of 
SAFDGS. 



RESEARCH DESIGN AND METHODS 

The SAFDGS family member recruitment and data collection procedures 
were reported previously (23). In brief, probands were Mexican Americans 
with T2DM with low income, and all first-, second-, and third-degree relatives 
of probands were invited to participate in the study. A variety of metabolic, 
hemodynamic, anthropometric, and demographic variables were collected 
from about 700 individuals drawn from 39 large Mexican American families. 
Diabetes status was defined by the 1999 criteria of the World Health Orga- 
nization (i.e., fasting glucose levels >126 mg/dL and/or 2-h glucose levels 
>200 mg/dL). 

Of the individuals examined at SAFDGS, GFR data were available for 
453 individuals distributed across 29 families. GFR values were estimated as 
described previously (5) by the two recently recalculated GFR prediction 
equations (24): the recalculated simplified four-variable Modification of Diet 
in Renal Disease (GFR-MDRD) study equation (25) and the Cockcroft-Gault 
(GFR-CG) equation (26) acyusted for body surface area and corrected for the 
bias in the MDRD data. The quantitative trait values of triglycerides (TGL) and 
albumin-to-creatinine ratio were log transformed and used in the association 
analyses since their raw data were nonnormally distributed. The Institutional 
Review Board of the University of Texas Health Science Center at San Antonio 
approved all procedures, and all subjects gave informed consent. 
Molecular variants identification and genotyping. Exons and 2-kb putative 
promoter region of IRSl gene were screened for DNA sequence variants in 32 
individuals, who contributed positively to the linkage, by direct sequencing of 
PGR products amplified from genomic DNA isolated from peripheral blood 
leukocytes. Oligonucleotide primers used for PGR and sequencing are avail- 
able upon request. Sequencing was performed using ABI Prism Big Dye Ter- 
minator Cycle Sequencing kit (Applied Biosystems [ABI]) and a capillary 
sequencer (Model 3730x1; ABI). Genotyping of all the SNPs was performed 
using the TaqMan assay (ABI), which was carried out on a GeneAmp PGR 
system 9700 (ABI), and fiuorescent signals were detected on an ABI PRISM 
7700 sequence detector (ABI). To assure accuracy of the genotyping, coded 
blind replicate samples from 50 subjects were included in each assay. 
Statistical genetic analysis. The genotypic data were checked for Mende- 
lian consistency using the program SimWalk2 (27). All polymorphisms were 
tested for Hardy-Weinberg Equilibrium. Allele frequencies were estimated 
using maximum likelihood techniques, which account for pedigree structure. 
Linkage disequilibrium (LD) between SNPs was estimated using the val- 
ues. We performed association analysis in our complex pedigree-based data 
using the measured genotype approach (MGA) within the variance compo- 
nents (VC) analytical framework (28,29). The VC-based approach accounts 
for the nonindependence among family members. In this analytical tech- 
nique, VCs are modeled as random effects (e.g., additive genetic effects and 
random environmental effects), whereas the effects of measured covariates 
such as age and sex are modeled as fixed effects on the trait mean. The 
marker genotypes were incorporated in the mean effects model as a mea- 
sured covariate, assuming additivity of allelic effects (28-30). The effect of 
this measured genotype (i.e., association parameter) together with other 
covariate effects and VCs were estimated by maximum likelihood tech- 
niques. For association analysis between genotypes and GFR, age and sex 
terms in addition to other covariates were used in consideration of both 
biological and/or statistical significance. For GFR-CG analysis, age, age^, sex, 
diabetes, duration of diabetes, systolic blood pressure (SBP), and antihy- 
pertensive medications were used as covariates. For GFR-MDRD, age, sex, 
BMI, diabetes, duration of diabetes, SBP, and antihypertensive medications 
were used as covariates. 

The hypothesis of no association is tested by comparing the likelihood of 
a model in which the effect of the measured genotype is estimated with a model 
where the effect of the measured genotype was fixed at zero. Twice the dif- 
ference in the log-likelihoods of these models yields a test statistic that is 
asymptotically distributed, approximating a distribution with one degree of 
freedom. Before MGA was performed, the quantitative transmission disequi- 
librium test (31) was used to examine hidden population stratification. Given 
the LD patterns among the examined 13 SNPs, we obtained a P value after 
adjusting for multiple testing using an effective number of SNPs calculated 
following the method of Li and Ji (32). For a statistical threshold of <0.05, the 
required experiment-wide significance threshold, using the effective number 
of SNPs (which was 7.66) to adjust for multiple testing was 0.00667. Because 
the issue of correction for multiple testing is not straight forward with regard 
to multiple correlated traits, there were no additional attempts to correct for 
multiple testing because the corresponding models are not independent. All 



statistical techniques described above were implemented in the computer 
program SOLAR. 

The Bayesian quantitative trait nucleotide (BQTN) analytical technique was 
used to analyze the IRSl SNP and GFR data further, with the use of the 
program SOLAR. This technique has been detailed elsewhere (33). Given 
complete SNP data for a gene, this statistical technique can be used to identify 
the sequence variants that are either potentially functional or that exhibit the 
highest disequilibria with such potential functional sites. The BQTN model is 
a simple extension of the classical variance component model, which aims to 
disentangle the genetic architecture of a quantitative trait. Because a candi- 
date gene may contain a number of SNPs that could generate several possible 
competing models of quantitative trait nucleotide (QTN) action, we used 
a Bayesian model selection/model averaging approach to analyze the SNP data 
simultaneously to estimate the probability that each SNP has a direct effect on 
the phenotype (33,34). 

To determine whether the associated SNP(s) found by the QTN analysis can 
account for our reported linkage of GFR at chromosome 2q35-37, we combined 
the QTN analysis with our identity-by-descent-based variance component 
linkage analysis. If a variant, or set of variants, in the IRSl gene contributes to 
the observed linkage signal, linkage analysis conditional on a fixed effect, 
MGA of the polymorphism will yield a LOD score near zero. Alternatively, if 
the associated polymorphism is in less than complete LD with the true func- 
tional site, linkage analysis will yield a reduced, but nonzero LOD score. This 
method has been described in detail elsewhere (29). 

Western blot analysis. The plasmids carrying IRSl gene with glycine at 
codon 972 (wild type) and arginine at codon 972 (mutant) were provided by 
Dr. Sesti (35). Plasmids were transfected into human kidney mesangial cells 
using the FuGene HD transfection reagent as recommended (Roche Diag- 
nostics). Cells (90-95% confiuence) were washed twice with PBS, and media 
were replaced with 0.45 mL of OPTI-MEM I (Invitrogen). Precomplex of the 
plasmid DNA (0.5 jjig/well) with FuGene (0.5 jjuL/well) in Opti-MEM was mixed 
(50 jjiL/well) and incubated at room temperature for 5 min. DNA and FuGene 
complexes were added to each well and incubated at 37°C with 5% CO2. Cells 
were then incubated for 6 h followed by addition of 0.5 mL of fresh media with 
34% serum to a final concentration of 17%. After 48 h of transfection, cells were 
made quiescent by serum deprivation for 24 h, followed by treatment with 
insulin (100 |jlM) for 0, 5, and 10 min. Cells were harvested and lysed with 
radioimmunoprecipitation assay buffer. Equal amounts of total cellular pro- 
teins were fractioned by SDS-polyacrylamide gel electrophoresis, transferred 
to nitrocellulose membranes, and probed with the following antibodies: anti- 
IRSl (Upstate), antiphospho-IRSl (Tyr 941) (Santa Cruz), anti-Akt, antiphospho- 
Akt (Thr 308) (Cell Signaling), and antiactin (Sigma). Bound antibodies were 
visualized using SuperSignal West Pico kit (Thermo Scientific). Densitometric 
analysis was performed using the National Institutes of Health ImageJ program. 

RESULTS 

The clinical characteristics of the study participants are 
shown in Table 1. Of the genotyped individuals, the avail- 
able phenotypic data for each phenotype varied from 610 
subjects for total cholesterol to 670 subjects for age. In the 
total sample, 29, 28, and 14% of the subjects had hyper- 
tension, T2DM, and albuminuria, respectively (Table 1). 
However, the GFR data were available for only 453 sub- 
jects. Of the individuals with GFR data, 29% had T2DM. 

IRSl (NM_005544) is composed of two exons, and the 
full-length protein is encoded by exon-1 (Fig. 1). To iden- 
tify DNA polymorphisms in IRSl, we sequenced both 
exons and 2-kb putative promoter region of IRSl in 32 
subjects who contributed positively to the linkage of GFR. 
In total, we identified 11 SNPs including seven in the 
coding region, and two each in the 3' UTR and putative 
promoter regions (Fig. 1). Of the seven SNPs identified 
in the coding region, five are synonymous (Asp90Asp, 
Gly234Gly, His400His, Ala804Ala, Pro893Pro) and two are 
nonsynonymous SNPs (Pro512Ala, Gly(972)Arg). In addi- 
tion, we chose eight tagging SNPs (SNPs-10-16, 19) from 
HapMap database for the analysis to comprehensively 
cover ~59-kb-long region of intron 1 of IRSl (Fig. 1). All 
the 19 SNPs (11 identified by resequencing and 8 from the 
public database) shown in Fig. 1 (SNP-1 to SNP-19) were 
genotyped in the entire data set (N = 670; 39 families) by 
TaqMan assay. 
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TABLE 1 

Clinical characteristics of the genotyped SAFDGS participants 
used for the current study^ 



Phenotypes Mean ± SD or % 





61 




TH/TiP Hijil^ptp<5 


28 




T-rvTiPrtpn <s;i on 


29 




AlVjiiminnria 


14 






44.8 ± 


16.2 


SBP rmmH^") 


128.6 ± 


18.2 


DBF (mmHg) 


71.3 ± 


9.6 


BMI (kg/m^) 


30.9 ± 


7.0 


Waist circumference (mm) 


1,001.0 ± 


170.6 


Total cholesterol (mg/dL) 


194.0 ± 


37.9 


HDL-C (mg/dL) 


46.0 ± 


12.0 


Triglycerides (mg/dL) 


152.7 ± 


100.3 


Log transformed triglycerides 


4.9 ± 


0.6 


GFR-CG (mL/min/1.73 m^) 


90.0 ± 


26.9 


GFR-MDRD (mL/min/L73 m^) 


88.4 ± 


17.9 


Albumin (mg/dL)-to-creatinine (mg/dL) ratio 


0.1 ± 


1.5 


Log transformed albumin-to-creatinine ratio 


2.4 ± 


0.8 



^Sample size varies from 610 (cholesterol) to 670 (age) except for 
GFR-CG and GFR-MDRD. The sample size for GFR-CG and GFR- 
MDRD = 453. 



Based on the genotypic data, six were excluded from 
further analysis because one (SNP-7) failed to be in Hardy- 
Weinberg Equilibrium and the minor allele frequencies of 
five SNPs (SNPs-5, 6, 8, 17, and 18) were less than 0.5% 
(Fig. 1). We estimated the pairwise LD (r^) between all 
the 13 SNPs and found that the pairwise LD ranged from 
0 to 1.0 (Fig. 2). As can be seen from Fig. 2, the relatively 
high pairwise LDs (r^ > 0.8) were found for the following 
SNPs pairs: rs2234931-rsl801278 (r^ = 1.0), rs2288586- 
rsl0498212 = 0.97), rs4675093-rs 10205923 (r^ - 0.92), 



rsl2052364-rsl0181778 (r^ = 0.88), and rsl2052364-rsl 1690232 
{f = 0.84). 

Association analysis between genotypic and GFR data 
were performed using MGA. The (B-coefficients and P val- 
ues of the covariates used are given in the Supplementary 
Data. The location, allele frequencies, and association 
analyses of 13 SNPs are summarized in Table 2. Of the 13 
SNPs examined for association, the G to A substitution at 
codon 972 changing the amino acids from Glycine to Ar- 
ginine (Gly(972)Arg) with minor allele frequencies of 
about 4% was found to be significantly associated with 
GFR-CG (P = 0.0006). Given the sample size of 431 for this 
analysis (Table 2 and Table 3), genetic effects that account 
for 2.9% of the total phenotypic variation in GFR-CG could 
be detected with 80% power with an experiment-wide 
P value < 0.05 (i.e., effective number of SNPs = 7.66). In 
this study, the Gly(972)Arg polymorphism explains 3.3% of 
the phenotypic variation in GFR-CG. The Gly(234)Gly 
(rs2234931), which was in complete LD (r^ = 1; Fig. 2) with 
Gly(972)Arg variant (rsl801278) also exhibited similar as- 
sociation with GFR. The Gly(972)Arg and Gly(234)Gly 
variants were also found to be significantly associated with 
GFR-MDRD (P = 0.0046) (Table 2). These associations 
continue to be significant after adjusting for multiple 
testing. It should be noted that no individuals were ho- 
mozygous for the rare allele in either of these SNPs. 

Subsequently, we conducted BQTN analysis with eight 
SNPs, after excluding five SNPs (SNP-4-rs2234931; SNP- 
12-rsl0498212, SNP-15-rsl0205923, SNP-16-rsl0181778, 
SNP-19-rsl 1690232) that exhibited redundancy owing to 
relatively high pairwise LD (r^ > 0.80) (Fig. 1 and Fig. 2). 
Of the eight SNPs, SNPs 1-3, 9-11, 13, and 14 resulted in 
the examination of 256 plausible additive gene action 
models; the Gly(972)Arg variant provided strongest evi- 
dence for the observed association with GFR-CG, with 
estimated posterior probability of >0.99 (data not shown). 
Such a high posterior probability of effect is suggestive of 
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FIG. 1. Schematic diagram of human (NM_005544) gene structure on 2q36 and the location of the SNPs identified and genotyped in SAFDGS. 
The exons of the IRSl are represented by solid box and the intron by a thin line. SNPs with the base change and their reference sequence (rs) 
numbers are indicated by the vertical arrow. 
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FIG. 2. LD between SNP pairs within the IRSl gene. SNPs are labeled on the ^-axis, and the locations (bp) within the gene are shown on the jcr-axis. 
Pairwise LD is estimated using values and depicted in the figure by the color intensity of the shaded box, as shown in the legend. The diagonal 
represents a comparison of each SNP against itself (i.e., = 1.0). (A high-quality color representation of this figure is available in the online 
issue.) 



a direct functional effect for this variant. As revealed by 
the MGA, the mean GFR values by genotypic categories 
differed significantly (Table 3). The carriers of Arg972 
variant had significantly decreased GFR-CG and GFR- 
MDRD values. 

To determine whether Gly(972)Arg could account for 
the initial linkage signal, we performed linkage analysis on 
chromosome 2 conditional on the measured genotype ef- 
fect. Because the sample size available (A^ = 431) for the 
Gly(972)Arg genotypic data was slightly smaller than the 
sample size (A^ = 453) used for the original linkage study 
(5), we repeated the linkage analysis of GFR-CG with 
genotype-by-diabetes interaction influences. As was done in 
our original linkage report (5), the LOD score obtained 
from the linkage analysis was adjusted with 2df to Idf, 
denoted as LODc score. More details on the adjustment of 
2df to Idf are provided in the Supplementary Data. The 
LODc obtained in the current study was 3.1 (A^ = 431) 
compared with the original report of 3.3 (A^ = 453). By 
inclusion of the Gly(972)Arg polymorphism in the model, 
the LODc score went down from 3.1 to 2.3; thus, about 26% 
of the evidence for linkage of GFR at 2q was explained by 



this variant (Fig. 3). Taken together, these data suggest that 
the Gly(972)Arg contribute to the decrease in GFR and that 
Gly(972)Arg may be one of the several potential functional 
sites regulating GFR or it is in LD with the true functional 
variant(s) in IRSl or genes flanking IRSl. 

To further determine the functional mechanism by 
which Arg972 of variant of IRSl may contribute to varia- 
tion in GFR, we transfected the human renal mesangial 
cells with the expression vector pcDNA3 carrying wild- 
type (Gly972) and mutant (Arg972) IRSl (35) and studied 
their effect on the insulin-induced tyrosine (941) phos- 
phorylation of IRSl and threonine (308) phosphorylation 
of Akt kinase, a critical downstream effector of PI 
3-kinase. Transfected cells were treated with 100 nmol/L 
insulin for 0, 5, and 10 min and lysed. Equal amount of total 
proteins were subjected to SDS-PAGE and probed with the 
following antibodies: anti-IRSl, antiphospho-IRSl, anti- 
Akt, antiphospho-Akt, and actin. In cells expressing Arg972 
variant of IRSl, insulin-stimulated phosphorylation of IRSl 
was significantly decreased by 57.3% at 5 min and by 65.8% 
at 10 min of insulin treatment compared with cells ex- 
pressing Gly972 variant of IRSl (Fig. 4A). Compared with 
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TABLE 2 

Association analysis between the IRSl polymorphisms and GFR and T2DM-related traits 



Traits^ and P values 



SNPs# 


OTVTT> TT\^ 

bJNP_lDs 


Position^ 


Location 


Alleles {%) 


HDL-C^ 


IGL 


bBP 


DBP^ 


(Jth K-CG 


(Jrh K-MDKD 


SNP-1 




-2990 


Promoter 


C(97.4)/T(2.6) 


0.0444 


0.6891 


0.1639 


0.1889 


0.5684 


0.5451 


SNP-2 


rs6436635 


-2221 


Promoter 


C(89.9)/T(10.1) 


0.2535 


0.0328 


0.0007 


0.0023 


0.7696 


0.7959 


SNP-3 


rs3731597 


270 


AspC90)Asp 


C(96.2)/T(3.8) 


0.0386 


0.8630 


0.4236 


0.6214 


0.6239 


0.8650 


SNP-4 


rs2234931 


702 


Gly(234)Gly 


G(96.5)/A(3.5) 


0.1555 


0.0033 


0.5016 


0.9860 


0.0006 


0.0044 


SNP-9 


rsl801278 


2911 


Gly(972)Arg 


G(96.5)/A(3.5) 


0.1494 


0.0033 


0.5299 


0.9719 


0.0006 


0.0046 


SNP-10 


rs3769647 


3872 


Intron-1 


C(84.4)/G(15.6) 


0.5946 


0.2234 


0.1537 


0.4142 


0.5016 


0.8834 


SNP-1 1 


rs2288586 


8065 


Intron-1 


C(81.6)/G(18.4) 


0.4114 


0.3943 


0.2982 


0.5108 


0.4556 


1.0000 










T(81.2)/A(18.8) 














SNP-13 


rs4675093 


25692 


Intron-1 


T(68.7)/A(31.3) 


0.0545 


0.0911 


0.0025 


0.0725 


0.5850 


0.6354 


SNP-14 


rsl2052364 


28247 


Intron-1 


T(83.2)/C(16.8) 


0.5167 


0.2311 


0.1673 


0.7027 


0.6893 


0.8596 1 


SNP-15 


rsl0205923 


48898 


Intron-1 


C(69.8)/T(30.2) 


0.0614 


0.2053 


0.0085 


0.1176 


0.5958 


0.6323 


SNP-16 


rsl0181778 


49218 


Intron-1 


T(82.4)/A(17.6) 


0.3172 


0.5441 


0.2514 


0.5857 


0.8832 


0.9968 1 


SNP-19 


rsl 1690232 


65614 


3' UTR 


A(83.9)/G(16.1) 


0.3639 


0.2971 


0.0790 


0.3825 


0.8143 


0.4537 



#See Fig. 1 for the location of each variant. ^Position in base pairs relative to the A of the start codon (ATG), which is +1. ^The availability of 
genotypic data together with the requirement for adjustment of covariate influences for a given trait resulted in different sample size for each 
of the traits reported in the table. However, the allele frequencies reported here correspond to the total available genotypic data. *^Adjusted for 
age and sex terms, diabetes, and lipid medication. ^Adjusted for age and sex terms, diabetes, and antihypertensive treatment. ^Adjusted for 
age and sex terms, diabetes, duration of diabetes, SEP, and antihypertensive treatment. ^Adjusted for age and sex terms, diabetes, duration of 
diabetes, SBP, BMI, and antihypertensive treatment. 



cells expressing Gly972 variant of IRS, insulin-induced 
Thr308-Akt phosphorylation was also significantly re- 
duced in cells expressing Arg972 variant of IRSl by 
75.9% at 5 min and by 88% at 10 min of insulin treatment 
(Fig. 4B). 

Association analysis was also extended to available 
phenotypic data including SBP, diastolic blood pressure 
(DBP), TGL, HDI^choIesteroI (HDL-C), albumin-to-creatinine 
ratio, and T2DM. It is noteworthy that the two variants [Gly 
(972)Arg and Gly(234)Gly] that exhibited significant asso- 
ciations with the GFR measures were also found to be 
significantly associated with TGL (P = 0.0033) (Table 2). In 
addition, a promoter variant of IRSl (SNP-2; rs6436635) 
with minor allele frequency of 10.1% was significantly as- 
sociated with TGL (P = 0.0328), SBP (P = 0.0007), and DBP 
(P = 0.0023) (Table 2). In addition, a promoter variant 
(SNP-1) with a minor allele frequency of 2.6% and a syn- 
onymous variant (SNP-3; AspQOAsp) with minor allele 
frequency of 3.8% also showed statistically significant asso- 
ciation with HDI^C (P = 0.0444 andP = 0.0386, respectively). 
A pair of intronic region SNP (rs4675093-rsl0205923), which 
were in strong LD (r^ = 0.92), also exhibited significant as- 
sociation with SBP (P = 0.0025 and 0.0085, respectively). 
These SNPs also showed suggestive association with HDLrC 
(Table 2). 



DISCUSSION 

There is a large interindividual variation in GFR in patients 
with diabetes and in the general population. It is known 
that variation in GFR is influenced by environmental and 
genetic factors as well as their interactive influences (2-11). 
In an effort to identify and characterize GFR susceptibility 
gene(s), we localized a genetic region on chromosome 
2q35-37 that exhibited linkage to GFR (5). As we report in 
this study, this locus, which does not confer susceptibihty to 
T2DM (5), exerts a differential influence on GFR in diabetic 
and nondiabetic individuals. Of the positional candidate 
genes, IRSl is a strong susceptibility candidate gene to influ- 
ence GFR since the IRSl protein plays a pivotal role in 
modulating tissue response to insulin including in the kid- 
ney. A structural defect in this protein as a result of genetic 
variation may alter kidney function. 

To date, several DNA sequence variants have been 
identifled in human IRSl. Of those, the main focus has 
been on GIy(972)Arg variant, since Almind et al. (36) first 
identified this polymorphism to contribute to T2DM in 
a Danish cohort. Subsequent case control studies in dif- 
ferent ethnic groups have examined the association be- 
tween the Gly(972)Arg variant of IRSl and T2DM and 
related traits. Although the association results have been 
controversial, a meta-analysis combining 27 studies, which 



TABLE 3 

Mean values of GFR by Gly(972)Arg genotype categories 







Genotype categories# 






Traits 


Gly/Gly 


N 


Gly/Arg 


N 


P values 


GFR-CG (mL/min/1.73 m^) (A^ = 431) 


84.77 ± 2.15 


396 


71.84 ± 4.15 


35 


0.0006 


Log transformed TGL (N = 588) 


5.02 ± 0.05 


550 


4.73 ± 0.05 


38 


0.0033 



#There were no homozygous individuals identifled for the minor allele (Arg972). The numbers of genotyped individuals available for the 
analyses of GFR-CG, GFR-MDRD, and log transformed TGL were 431, 429, and 588, respectively. The requirement for adjustment of covariate 
influences for a given trait resulted in slightly reduced sample size for each of the three traits. 
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FIG. 3. Linkage analysis conditional on the Gly(972)Arg polymorphism on chromosome 2q. 



included 3,408 cases and 5,419 control subjects, suggested 
that the carriers of the Gly(972)Arg variant have a 25% 
increased risk for developing T2DM (17). In addition, 
a marked ethnic difference in the distribution of Gly(972) 
Arg variant was also observed. The Arg972 variant of IRSl 
was found to be more common in Caucasians (20%) (37) 
compared with African Americans (11%) (38), Asians (4%) 
(39), and Mexican Americans (4%) (40). In this study, we 
observed a frequency of about 4% for Arg972, which is 
similar to the one reported in another Mexican American 
population (40). Similarly, Gly(972)Arg failed to show an 
association with T2DM as reported (40). 

Of the 13 SNPs examined for association with GFR, only 
the Gly(972)Arg variant of IRSl exhibited significant as- 
sociation with GFR estimated by two formulae. The BQTN 
analysis using information from eight nonredundant SNPs 
also demonstrated that the Gly(972)Arg variant could in- 
fluence phenotypic variation in the levels of GFR with 
estimated posterior probability of >0.99. This high poste- 
rior probability effect reflects the potential functional im- 
pact of this variant on GFR. It is interesting that the 
carriers of Arg972 had significantly decreased GFR values. 
The subsequent conditional linkage/QTN analysis revealed 
that the Gly(972)Arg variant explained —26% of the linkage 
signal for GFR on 2q. Taken together, our data demon- 
strate that the Gly(972)Arg variant significantly contributes 
to decrease in GFR. 

The functional mechanism by which Arg972 variant 
contributes to decline in GFR remains speculative. Several 
studies have demonstrated that expression of the 972Arg 
variant of IRSl in cultured adipocytes, muscle, and pan- 
creatic b cells caused a defect not only in Tyr-941 phos- 
phorylation of IRSl but also in binding of the p85 regulatory 
subunit of PI3K to the IRSl, resulting in a decrease in 
IRSl-associated PI3K activity and subsequent activation 
of the Ser/Thr kinase Akt (19,35,41-44). Hribal et al. (20) 
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generated transgenic mice overexpressing Arg972 and 
confirmed a significant attenuation in insulin signaling path- 
ways downstream of IRSl in the liver, skeletal muscle, and 
adipose tissue. 

In the current study, we studied the effect of Arg972 on 
the insulin-stimulated tyrosine phosphorylation of IRSl 
and threonine phosphorylation of Akt in human renal 
mesangial cells. We show that expression of Arg972 in 
mesangial cells significantly reduces insulin-induced phos- 
phorylation of IRSl and Akt kinase compared with wild- 
type controls. Impaired IR/IRSl signaling in mesangial cells 
provides a potential mechanism by which Arg972 variant 
of IRSl contributes to decreased GFR. Alterations in the 
contractile property of mesangial cells as a result of im- 
paired IR/IRSl signaling by Arg972 variant may lead to 
alterations in glomerular hemodynamics, thereby influenc- 
ing GFR. Because insulin has been shown to be a survival 
factor for mesangial cells (45), impaired insulin signaling by 
the Arg972 variant may lead to apoptosis of these cells, 
thereby contributing to structural changes in the glomerulus 
resulting in a decrease in surface area available for filtra- 
tion. Functional studies are necessary to prove the hy- 
pothesis that the Arg972 variant of IRSl not only impairs IR 
signaling to protect mesangial cells against apoptosis but 
also impairs contractile function of mesangial cells. Im- 
paired nitric oxide (NO) production in vascular endothelial 
cells by the Arg972 variant may also contribute to the de- 
creased GFR. Human umbilical vein endothelial cells iso- 
lated from carriers of the Arg972 variant of IRSl have 
impaired eNOS expression in response to insulin (46). Be- 
cause the impaired IR/IRSl signaling is a potential mecha- 
nism for the observed association, we also performed 
association analysis for the 13 SNPs using fasting glucose 
data (A^ = 541) after excluding T2DM individuals that are 
on antidiabetic medications. Although the Gly(972)Arg 
variant failed to show an association, the promoter variant 
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FIG. 4. Effect of Arg972 on the insulin-induced phosphorylation of IRSl 
and Akt in human mesangial cells. A: Tyrosine (941) phosphorylation of 
IRSl. Human mesangial cells were transfected with plasmids carrying 
wild type (WT; Gly972) and mutant (Mut; Arg972) IRSl and treated 
with 100 nM insulin for 0, 5, and 10 min. Equal amount of total cellular 
proteins were separated on SDS-PAGE gel and probed with anti-IRSl 
and antiphospho-IRSl. B: Threonine (308) phosphorylation of Akt ki- 
nase. Human mesangial cells were transfected with plasmids carrying 
wild-type (Gly972) and mutant (Arg972) IRSl and treated with 100 nM 
insulin for 0, 5, and 10 min. Equal amount of total cellular proteins were 
immunoprecipitated with anti-Akt kinase and separated on SDS-PAGE 
gel and probed with anti-Akt and antiphospho-Akt kinase. Represen- 
tative images of three independent experiments are shown. 



perform (e.g., inulin, iothalamate clearance) in large family 
studies, the estimated GFR using the abbreviated MDRD 
and GC formulae has been proposed as the best validated 
means for transforming serum creatinine measurements 
into GFR (25,51) and used in several linkage studies (2-11). 
Another limitation is that because our resequencing effort 
was based on only 32 individuals, it is possible that we 
might have missed rare variants specific to this population. 
However, given the recent emphasis on the role of rare 
variants in complex disease causation compared with the 
common variants examined using the genome-wide asso- 
ciation studies, it is very interesting to note that the Gly 
(972)Arg associated with decline in GFR values has a fre- 
quency of 4.0%. We tried to replicate the association between 
Gly(972)Arg variant of IRSl and GFR in two independent 
family studies: the San Antonio Family Heart Study (Mexican 
Americans, N = 844) and the Joslin Study on the Genetics of 
Type 2 Diabetes (Caucasians, N = 825). Unfortunately, we 
failed to replicate the findings. 

In conclusion, we provide evidence for the first time that 
Arg972 variant of IRSl is associated with variation in GFR; 
carriers of Arg972 had significantly decreased GFR values 
and Arg972 variant accounts for 26% of the linkage of GFR 
on 2q37. It is noteworthy that we show that expression of 
Arg972 in human mesangial cells reduces the insulin- 
stimulated phosphorylation of IRSl and Akt kinase. To- 
gether, it is possible to speculate that the associated variant 
is one of multiple potential functional variants infiuenc- 
ing variation in GFR in the IRSl gene or is in LD with 
a potential functional variant(s) in the genes fianking 
IRSl. There is evidence for long-range LD in the human 
genome (52). However, in the absence of knowledge on 
such long-range LD of the SNP of interest in this study, it 
is important to note that genetic variants located 500 kb 
upstream of IRSl were recently shown to be strongly as- 
sociated with impaired metabolic profiles including insulin 
resistance and T2DM. Furthermore, these variants also im- 
paired the IR/IRSl signaling (21,22). Therefore, a thorough 
screening of the critical region on 2q37 may identify ad- 
ditional variants that could potentially relate to the GFR 
linkage finding in this study. Identification of susceptibility 
genes that contribute to the decline in GFR may facilitate 
prediction, prevention, and development of improved 
treatments of this devastating complication of diabetes. 



(rs6436635) was found to be nominally associated (P = 
0.035) with glucose levels. 

The present data also show an association of GIy(972)Arg 
variant with TGL in our cohort, and the carriers of Arg972 
had decreased TGL values. Our finding of lower TGL con- 
centrations in carriers of Arg972 contradicts previous 
observations in other populations at high risk to develop 
atherosclerotic cardiovascular disease (47,48). However, 
our findings are in agreement with others who reported 
that patients with T2DM (49) and patients with insulin 
resistance (50) have low concentrations of TGL. The 
functional significance of promoter variants, Asp90Asp, 
and a pair of intronic variants (rs4675093 and rs 10205923) 
that exhibit an association with cardiovascular risk factors 
such as HDL-C, TGL, SBP, and DBF in this study remains 
to be determined. 

Despite several strengths, our study has potential limi- 
tations. To our knowledge, the two estimates of GFR 
(GFR-CG and GFR-MDRD) have not been validated in the 
Mexican American population. Whereas the time-consuming 
and expensive direct measurement of GFR is difficult to 
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